Lagrangian theory of structure formation in relativistic cosmology II: 
average properties of a generic evolution model 
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Kinematical and dynamical properties of a generic inhomogeneous cosmological model, spatially 
averaged with respect to free-falling (generalized fundamental) observers, are investigated for the 
matter model 'irrotational dust'. Paraphrasing a previous Newtonian investigation, we present 
a relativistic generalization of a backreaction model based on volume-averaging the 'Relativistic 
Zel'dovich Approximation'. In this model we investigate the effect of 'kinematical backreaction' on 
the evolution of cosmological parameters as they are defined in an averaged inhomogenous cosmology, 
and we show that the backreaction model interpolates between orthogonal symmetry properties 
by covering subcases of the plane-symmetric solution, the Lemaltre-Tolman-Bondi solution and 
the Szekeres solution. We so obtain a powerful model that lays the foundations for quantitatively 
addressing curvature inhomogeneities as they would be interpreted as 'Dark Energy' or 'Dark Matter' 
in a quasi-Newtonian cosmology. The present model, having a limited architecture due to an 
assumed FLRW background, is nevertheless capable of replacing 1/4 of the needed amount for 
'Dark Energy' on domains of 200 Mpc in diameter for typical (one-sigma) fluctuations in a CDM 
initial power spectrum. However, the model is far from explaining 'Dark Energy' on larger scales 
(spatially), where a 6% effect on 400 Mpc domains is identified that can be traced back to an on 
average negative intrinsic curvature today. One drawback of the quantitative results presented is the 
fact that the epoch when backreaction is effective on large scales and leads to volume acceleration 
lies in the future. We discuss this issue in relation to the initial spectrum, the 'Dark Matter' 
problem, the coincidence problem, and the fact that large-scale 'Dark Energy' is an effect on the 
past light-cone (not spatial), and we pinpoint key-elements of future research. 
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I. INTRODUCTION 

In a previous work [LI] (henceforth cited as [LI]) we 
developed a Lagrangian theory for the evolution of struc- 
ture in irrotational dust continua in the framework of 
general relativity. This theory is characterized by em- 
ploying as a single dynamical variable the set of spatial 
coframes rf"^ , in the comoving-synchronous metric form: 



9ij 



with the spatial metric coefficients gij — Sab V^'iV^j, where 
are Gaussian normal (Lagrangian) coordinates that 
are constant along flow lines (here geodesies); i,j,k = 
1,2,3 denote coordinate indices, while a,b,c = 1,2,3 
are just counters. Having written Einstein's equations 
in terms of the variable 77°, we are entitled to set up 
a perturbation scheme that only perturbs this variable, 
while other variables like the density, the metric or the 
curvature, are expressed as functionals of this perturba- 
tion. Such a point of view extrapolates these latter phys- 
ical quantities beyond a strict perturbative expansion. In 
this spirit we defined and evaluated the first-order per- 
turbation scheme on a FLRW background cosmology as a 
clear-cut definition of what has been known as the "Rel- 



ativistic Zel'dovich Approximation" (henceforth abbrevi- 
ated by RZA) in previous work (see [LI] and references on 
previous work therein, in particular [37, 47, 49, 59, 60]). 

In the present work we will model the fluctuations us- 
ing this approximation and combine the resulting back- 
reaction model with exact average properties of inhomo- 
geneous dust cosmologies. (This framework is reviewed 
in [15] and [16] and will be briefly recalled below.) We 
will so paraphrase a detailed investigation in Newtonian 
cosmology [BKS] (henceforth cited as [BKS]), where we 
presented the results of an in-depth quantitative study 
of the kinematical backreaction effect in Newtonian cos- 
mology [10]. In order to evaluate the influence of inho- 
mogeneities on effective cosmological parameters we em- 
ployed in [BKS] exact solutions with planar and spherical 
symmetry, as well as perturbative solutions for generic 
initial conditions in the Eulerian and Lagrangian first- 
order perturbation schemes. The most general, non- 
perturbative model studied was based on the first-order 
Lagrangian deformation that was used as input into the 
kinematical backreaction functional (fluctuation model) , 
but the effective cosmology (the average) was then eval- 
uated within the framework of the exact averaged equa- 
tions. We will here provide the corresponding model in 
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general relativity. 

In [BKS] our initial data setting was conservative, i.e. 
we started with a standard Cold Dark Matter power 
spectrum without a cosmological constant at the begin- 
ning of the matter-dominated era, maintaining the Early 
Universe description as in the standard model. Perhaps 
the most surprising outcome of this study was that, al- 
though the kinematical backreaction term Q-p itself was 
quantitatively negligible on large expanding domains, the 
other cosmological parameters could still experience sig- 
nificant changes, e.g. for the deceleration parameter we 
found typical fluctuations up to 30 percent on the scale 
of 200 Mpc (with h = 0.5). In this work we will give the 
results for two background models: (i) a CDM spectrum 
with h = 0.5 as in [BKS] to allow for comparison, and (ii) 
a ACDM spectrum with h = 0.7 to study the effect of the 
background. For the interpretation of the emerging cur- 
vature inhomogeneities as 'Dark Energy' we refer to the 
first model. Note that both homogeneous-isotropic mod- 
els are assumed to have Euclidean space sections. Recall 
that, in Newtonian cosmology, the kinematical backreac- 
tion vanishes by construction on the scale of an imposed 
torus architecture (i.e. the use of periodic boundary con- 
ditions) [10], a restriction that no longer applies to the 
general relativistic model. The present investigation thus 
offers a more general view on the backreaction model, no- 
tably by allowing for a nontrivial curvature evolution. 

We may approach the present work with the following 
guidelines: 

• The cosmic triangle [3] of the three parameters of 
the standard model (related to the homogeneous 
density, the constant curvature and the cosmolog- 
ical constant) fluctuates on a given spatial scale. 
The measure of deviations from the FLRW time- 
evolution of these parameters is a kinematical back- 
reaction functional Q-p that encodes averaged fluc- 
tuations in kinematical invariants (the rate of ex- 
pansion and the rate of shear). 

• The influence of backreaction may be two-fold: for 
dominating shear fluctuations the effective density 
decelerating the expansion is larger than the ac- 
tual matter density, thus mimicking the effect of 
a Kinematical Dark Matter, for dominating ex- 
pansion fluctuations backreaction acts accelerating, 
thus mimicking the effect of a (positive) cosmolog- 
ical constant, i.e. Kinematical Dark Energy. 

• In general relativity, backreaction entails an emerg- 
ing intrinsic curvature. Thus, as soon as we require, 
on some scale, that the average model is given by a 
standard FLRW model geometry, then the emerg- 
ing curvature term allows to interpret the above- 
mentioned kinematical effects as 'Dark Matter' or 
'Dark Energy', respectively. In Newtonian cosmol- 
ogy, the backreaction term should be interpreted as 
cosmic variance of velocity fluctuations on a Euclid- 
ian space; a general-relativistic model entails devi- 
ations from an Euclidian model geometry; we so are 



led to interpret the dark sources in an assumed Eu- 
clidian background as "curvature energies" of the 
actual average distribution. 

• A genuinely relativistic property is the coupling of 
the kinematical backreaction functional to the av- 
eraged scalar (3-Ricci) curvature. This coupling 
is absent in weakly perturbed FLRW spacetimes, 
e.g. in a Newtonian framework where a fluctuating 
scalar curvature is absent, or in a post-Newtonian 
framework and standard relativistic perturbation 
theory at first order, where the cosmological model 
is required to stay close to the constant-curvature 
FLRW spacetime, if this latter is taken as the back- 
ground spacetime. 

• If kinematical backreaction is positive, i.e. domi- 
nated by expansion fluctuations and by a negative 
averaged curvature, a small source term is capa- 
ble of driving the averaged system away from the 
homogeneous-isotropic FLRW background. We 
found, with the help of a detailed dynamical sys- 
tems analysis of scaling solutions [57] , generalizing 
a previous analysis [18], that the standard model 
is a repellor in this situation, and the physical 
background (the average) is no longer given by the 
FLRW background. We may say that the physical 
background lies in the 'Dark Energy sector'. A sim- 
ilar situation occurs in the case of shear-dominated 
evolution of small-scale domains with positive av- 
eraged curvature, where the average is driven to 
lie in the 'Dark Matter sector'. This insight sug- 
gests, that a perturbation theory has to be set up 
on a background that interacts with structure for- 
mation (a framework for such a perturbation theory 
has been suggested [58]). This remark implies that 
the present investigation is limited, since we employ 
perturbations on a FLRW background. The result- 
ing backreaction model, however, can be formally 
employed also in such a more general framework, 
which is the subject of forthcoming work. 

We proceed as follows. In Sect. II we briefly review the 
equations governing the averaged inhomogeneous dust 
cosmologies in the relativistic framework and provide dif- 
ferent representations of the equations that are used in 
this article. We then move in Sect. Ill to dynamical 
models that allow estimating the backreaction terms: we 
provide backreaction and curvature models for the La- 
grangian linear perturbation theory, restricted to a rel- 
ativistic form of Zel'dovich's approximation investigated 
in detail in [LI]. Sect. IV is devoted to consistency checks 
including the Newtonian limit. The plane, spherical and 
quasi-spherical Szekeres solutions serve as exact refer- 
ence models in Sect. V. We then illustrate the results in 
Sect. VI in terms of the time evolution of cosmological 
parameters: we examine the evolution of the volume scale 
factor for typical initial conditions; then we explore the 
cosmological density parameters in the averaged model. 
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Here, we also quantify the emerging intrinsic curvature. 
We discuss the resuhs and conclude in Sect. VII. 



II. AVERAGED DUST MODEL IN GENERAL 
RELATIVITY 

A. Averaging 

In this section we will briefly set up the averaging 
framework developed in [12-14]. 

Given a foliation of spacetime into flow-orthogonal hy- 
persurfaces, we investigate fluid averaging for the matter 
model 'irrotational dust'; the geometry of the dust con- 
tinuum is described by the spatial metric coefficients gij 
in the comoving and synchronous line-element: 
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dX'dX^ 



(1) 



where X' are the already introduced (Lagrangian) co- 
ordinates that are constant along flow lines (spacetime 
geodesies of freely falling observers). Proper time- 
derivative along the 4— velocity u^, that is here equal 
to the coordinate time-derivative, will be denoted by 
df '■= u^d^ and sometimes by an overdot (where Greek 
indices label spacetime and Latin ones spatial coordi- 
nates). 

We look at the scalar parts of Einstein's equations, av- 
eraged over a compact spatial domain 2? at time t that 
evolved out of the initial domain at time ti , conserv- 
ing the material restmass inside the domain. The volume 
scale factor ou, depending on content, shape and posi- 
tion of the domain V, is defined via the domain's volume 
Vvit) = \V\, and the initial volume Vd, = V-DiU) = \Vi\: 



1/3 



(2) 



where the volume of the domain is given by: 



Vv{t) := f d^l 
Jv 



(3) 

Jd^X, and 



with the Riemannian volume element d/ig 
the local volume deformation J :— det {gij ) / det (G^ j ) , 

Gij = 9ij{ti)- 

For domains V, which keep the total material mass 
Mx> constant, as for Lagrangian domains with geodesic 
motion of their boundaries, the average density is simply 
given by: 



M- 



V 



Mr 



Mt 



(4) 



In this setting we can derive equations analogous to the 
averaged equations derived in Newtonian cosmology [10] 
(restricting attention to irrotational dust flows). Spatial 
averaging of a scalar field 'J is defined by: 



(vI'(t,X^))^:=-l / d^^, *(t,X^) 



(5) 



The key-property of inhomogeneity of the field ^' is re- 
vealed by the rule of non-commutativity [10, 12]: 



(©)i,Wi, 



dt {'fit, x'^))^ - (dt^it, x'))^ - (e*)^ - 

(6) 

where O denotes the trace of the expansion tensor. In the 
present spacetime setting the latter is, up to the sign, the 
extrinsic curvature tensor Qij = —Kij. Q describes the 
local rate of volume change of fluid elements along the 
vector field dt- 

It is formally convenient in the following calculations, 
but not necessary, to introduce the following formal av- 
erage: 



1 



(A), = — I d^X A 



(7) 



where Vp, /p. d^X = J{X\ti)d^X is the volume 

of the initial domain Vi for J(X',ti) = 1 (we will ex- 
plain below that this latter property of the local volume 
deformation holds for non-normalized coframes). This 
average coincides with the Riemannian volume average, 
if we consider fields at initial time, A = A{ti) (therefore 
we label it by I). Rewriting the Riemannian average 
in terms of this formal average operation, we obtain the 
useful formula: 



(A), 



{J)i 



'■v ■ 



(8) 



It will simplify the following calculations and is useful to 
discuss the Newtonian limit (see Section IV A) . 



Equations for the evolution of average 
characteristics 



1. Generalized Friedmann equations 

The spatially averaged equations for the volume scale 
factor ax), respecting restmass conservation within the 
domain T), read ]12]: 
averaged Raychaudhuri equation: 



3^+4^G 



Mr 



- A = 



averaged Hamilton constraint: 



av 



3 V-DiO: 
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Qv 
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(9) 



(10) 



where the total restmass A/pj, the averaged spatial 
3-Ricci scalar (TVjj, and the kinematical backreaction 
Q-p are domain-dependent and, except the mass, time- 
dependent functions. The kinematical backreaction 
source term itself is given by: 



Q 



V • = 



2(11) 



V 



V 



(©). 
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(11) 
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Here, I and II denote the principle scalar invariants (de- 
fined equivalently to Eq. (44) below) of the extrinsic cur- 
vature coefficients Kij. The second equality follows by 
introducing the decomposition of the extrinsic curvature 
into the kinematical variables, by means of the identity 

0V := — The kinematical variables are the rate of 
J J ... 

expansionQ := —K^^. and the shear a^: = — JiT' — -^OfJ* -. 



We also defined the rate of shear := 1/2 a'^^a''^. 
We appreciate a close correspondence of the GR equa- 
tions (9) and (10) with their Newtonian counterparts 
(see [BKS]). The first equation is formally identical to 
the Newtonian one, while the second delivers an addi- 
tional relation between the averaged scalar curvature and 
the kinematical backreaction term that has no Newto- 
nian analogue. This implies an important difference that 
becomes manifest by looking at the time-derivative of 
Eq. (10). The integrability condition that this time- 
derivative agrees with Eq. (9) is nontrivial in the inho- 
mogeneous GR context and reads: 

dtQv + GHvQv + dt (7^)p + 2Hv (7^)p = , (12) 

where H-p ■— a-p/a-p denotes the volume Hubble func- 
tional. 

Eq. (12) shows that averaged scalar curvature and the 
kinematical backreaction term are directly coupled, un- 
like in the Newtonian case, where the curvature is non- 
existent. For initially vanishing k in the standard FLRW 
cosmology, the scalar curvature remains zero. This is 
not the case in the inhomogeneous GR context, where 
kinematical backreaction produces averaged curvature in 
the course of structure formation, even for domains that 
are on average flat initially. (Note also that the sign of 
the averaged curvature may change during the evolution 
contrary to FLRW models.) To express the deviation 
of cosmic curvature from a constant-curvature (quasi- 
Friedmannian) behavior, we define the average peculiar- 
scalar curvature by 

Wv ■■= {n)^-6kvjal . (13) 

Integrating Eq. (12) with this definition inserted, we ob- 
tain: 

5i(f ^°^<''')- 5 + 

(14) 

with c-pj = + Q-Dj. Inserting it back into Eq. (10) 
we obtain the generalized Friedmann equation: 

,2 




SttG {q) 



V 



A 
3" 



(15) 



This equation is formally equivalent to its Newtonian 
counterpart [10]. It shows that, by eliminating the aver- 
aged curvature, the whole history of the averaged kine- 
matical fluctuations acts as a source of a generalized 



Friedmann equation. This equation was the starting 
point of our investigations in [BKS] . Since it is also valid 
in general relativity, we are in the position to translate 
many results from ]BKS] into the GR context. In partic- 
ular, our numerical codes can be accordingly applied. 

We now provide a compact form of the averaged equa- 
tions introduced above, as well as some derived quantities 
that we will analyze in this paper. 



2. Effective Friedmannian framework 

The above equations can formally be recast into stan- 
dard Friedmann equations for an effective perfect fluid 
energy momentum tensor with new effective sources ]13]: 
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(16) 
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(17) 



Eqs. (17) correspond to the equations (9), (10) and (12), 
respectively. 

This system of equations does not close unless we im- 
pose a model for the inhomogeneities. Note that, if the 
system would close, this would mean that we solved the 
scalar parts of the GR equations in general by reduc- 
ing them to a set of ordinary differential equations on 
arbitrary scales V. Closure assumptions have been stud- 
ied by prescribing a cosmic equation of state of the form 
p^ff = (3{gfg,axi) [14, 15[, or by prescribing the backre- 
action terms through scaling solutions, e.g. Qx> oc 
[18; 57[. In this paper we are going to explicitly model 
Qt) by a relativistic Lagrangian perturbation scheme. 



3. The cosmic quartet and derived parameters 

For later convenience we introduce the set of dimen- 
sionless average characteristics: 



ttG 

— i0)v 



3H. 



— - 



:= 



A 



Q 



V 



(18) 



We will, henceforth, call these characteristics 'parame- 
ters', but the reader should keep in mind that these are 
functionals on T). Expressed through these parameters 
the averaged Hamilton constraint (10) assumes the form 
of a cosmic quartet [11[: 



V 



1 



(19) 
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In this set, the averaged scalar curvature parameter and 
the kinematical backreaction parameter are directly ex- 
pressed through (7?.)23 and Qt>, respectively. 

In order to compare these parameters with the 'Fried- 
mannian curvature parameter' that had to be used in 
Newtonian cosmology [BKS] , and that is employed to in- 
terpret observational data, we can alternatively introduce 
the set of parameters 



2 TT2 



— 5 i-j^ J 4-yy 



ml 



being related to the previous parameters by 



(20) 



(21) 



Like the volume scale factor a-p and the volume Hubble 
functional H-p, we may introduce 'parameters' for higher 
derivatives of the volume scale factor, e.g. the volume 
deceleration functional 



1 



V 



2nl 



(22) 



(For higher derivatives such as the statefinders, see [15].) 
In this paper we denote all the parameters evaluated at 
the initial time by the index 2?i, and at the present time 
by the index Vq. 



kinematical backreaction term of Eq. (11) in terms of the 
spatial coframe coefficients . These latter are defined 
through the metric form coefficients 



(23) 



where Gab^^^^j = Gij ~ gtj (ti,X'^^ are the initial met- 
ric coefhcients. The coframes therefore contain the com- 
plete time evolution of the metric. We now express the 
expansion tensor, 

with the help of the coframes. Since the inverse met- 
ric is analogously decomposable into frames (that are 
inverse to the coframes £^^'77" = ; e^r]\ = S^), 



we can finally write (for details see [LI]): 

p)i — p i^-,a _ ^ , JM -a b c 
^ j ^ ^aV j — r,j ^abc£ V jV kV e 



(24) 



(25) 



with the local volume deformation (corresponding to the 
Jacobian in Newtonian theory) 



1 



J■.= deti^^1)^-eabce'''v''^VJV'k 



(26) 



This may then be used to calculate the first two invariants 
of J and to obtain the expression 



III. QUANTIFYING BACKREACTION 

In the last section we saw how the appearance of the 
backreaction term modifies the standard Friedmannian 
evolution. In the following we will quantify this depar- 
ture. Exact inhomogeneous solutions for estimating the 
amount of backreaction are only available for highly sym- 
metric models like for models with plane symmetry (Sub- 
sect. V A) and the spherically-symmetric solutions (Sub- 
sect. VB). In generic situations we have to rely on ap- 
proximations. As the Newtonian result [BKS[ suggests, 
we will obtain also here an approximation that interpo- 
lates in special cases between the above two subclasses 
of exact GR solutions with orthogonal kinematical prop- 
erties. However, the classes of solutions will be more 
strongly restricted than in the Newtonian case, because 
of the Hamilton constraint that joins the system of equa- 
tions in the GR context. 



A. Modeling the local deformation 

1. Representation of deformations through coframes 

To be able to use the 'Relativistic Zel'dovich Approx- 
imation' as formulated in [Ll[, we will first express the 




3 (J) 



, (27) 



for Qd that solely depends on the coframes ry'J and their 
time-derivatives. 



2. The 'Relativistic Zel'dovich Approximation' 

The equations and 'parameters' introduced in Sec- 
tion II can live without introducing a background space- 
time. The description is background-free and non- 
perturbative. Also expressing Qx> in terms of coframes 
as sketched in the previous section is still completely gen- 
eral. However, in order to provide a concrete model for 
the backreaction terms, we employ methods of pertur- 
bation theory. We will, however, only model the fluc- 
tuations by perturbation theory; by using the exact ex- 
pressions for the functionals in terms of coframes, we 
extrapolate the first-order perturbative solution. 

We now introduce perturbed coframes analogous to 
[Ll[. In perturbation schemes one usually defines a ref- 
erence frame through a known solution, e.g. a frame 
comoving with the Hubble flow (i.e. a FLRW solu- 
tion) that we now represent through three homogeneous- 
isotropic deformation one-forms (labelled by a, &, c... and 



6 



expressed in the local exact coordinate basis labelled by 
i,j,k...), 

tJh^ - r^H^ dX^ a(t)T,^"(ti) , r^^^ := a{t)d 'J , (28) 

where a(t) is a solution of Friedmann's differential equa- 
tion 



^2 ^87TGgH{t)+A k 



3 



it) 



(29) 



with H{t) := a/a, QH{t) ~ QHja?{t) and a{ti) = 1. For 
the full deformation one-forms we prescribe the super- 
position 



77'^-77^'^ + a(^)P^ 



(30) 



with inhomogeneous deformation one-forms 'P°'{t,X^). 
To first order, they can be restricted to the relativistic 
generalization of Zel'dovich's approximation (RZA) [LI]: 

P" = ^{t)P1 dX^ , (31) 

with Pi := Pi {ti,X''). The function ^ (t) is defined by 

m m - q{u)]/m) , (32) 

where the function q{t) is solution of the equation: 

«,) + 2|1#)+(3^|-a)„„.„^ (3.3) 
Thus, the function ^ satisfies: 

(34) 

Explicit solutions of this equation for different back- 
grounds may be found in [9], those including a cosmo- 
logical constant in [4]. 

Writing the RZA of Eq. (30) in the exact coordinate 
basis dX' yields: 

R^Ary'^, (t, X'^) := a (t) (iV^ + ^(t)/"',) , (35) 

where N"", "^^ry" {tuX'') and P'\ = P° (ti, X'^). 

Before we use this expression to determine the backre- 
action term, let us elaborate on a subtlety with the choice 
of initial conditions. The expression of the metric tensor 
in terms of nonintegrable coframes: 



gtj Gab ifyfj , 



(36) 



allows two different treatments of the initial displace- 
ments. One can either include them into Gab which 
means that the noncoordinate basis is orthogonal but 
not orthonormal. Or one can choose Gab to be orthonor- 
mal (being the standard assumption), i.e. Gab = ^at, 
but then one has to deal with the initial values of the 
coframes. To have both at a time, i.e. N°-i = and 
Gab = Sab, IS uot possiblc as this would mean that the 



RZA initial metric would be Euclidian, and this would 
disable any time-evolution of this metric as pointed out 
in [47] and [59]. As a nontrivial time-evolution of the 
metric is what we are interested in, there are only two 
options: 

01: if iV", := + P''^, with P", := P%{U,X^) 
and P"-i 7^ 0, then G can be restricted to 
(5, and the coframes read (i, := 

a (t) i^S"^ + P"., + £,{t)P%^ with the metric g^j := 

SabV^'.v). 

02: by appropriate coordinate transformations, one 
may set N°-i in Eq. (35) to 5"; the transforma- 



tion then sends — )• 



(U,X'^)P», 



~ i . .. ^ 

and all information about the initial geometrical 
inhoniogeneities is contained in G. The coframes 



become RZA^a := a (t) (^(5° + ^(t)^° j , and 

the metric := Gab 'il'i'n''j with Gab = 5cd N'^aN'^b, 
with iV"^ as defined above. 

In order to have a complete formal correspondence with 
the Newtonian model, we stick to the second option in 
what follows. (Notice that we have chosen the first option 
in [Ll[.) 



B. Kinematical backreaction and intrinsic 
curvature models 

After the definition of the perturbative setup explained 
in the previous section, we now turn to the concrete cal- 
culation of the backreaction and curvature models. As 
discussed in the previous section, we use non-normalized 
coframes (Option 2; the corresponding expressions with 
Option 1 will be listed in appendix B), 



(37) 



As pointed out by Chandrasekhar [22[ such a choice can 
lead to formally simpler expressions in some cases, and 
here we encounter such a case. Indeed, the expression for 
the kinematical backreaction term turns out to resemble 
more closely its Newtonian counterpart, and the formulae 
become more concise. 

The 'Relativistic Zel'dovich Approximation' is then de- 
fined as: 



RZA 



ri^,{t,X''):=a{t)i5'^^+m 



with: 



a(ti) = 1 . 



(38) 



By definition (see [Ll[), the RZA of any field is the evalu- 
ation of this field in terms of its functional dependence on 
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the RZA coframes, Eq. (38), without any further approx- 
imation or truncation. Therefore, the metrical distances 
are calculated exactly for the approximated deformation: 



RZA 



7>b 



+e'WGa&^",^''j} , (39) 

and so 

'^^^.,(ti)=:G,, . (40) 
The local volume deformation, Eq. (26), then becomes: 



RZA 



(41) 



(implying that for the approximated deformation we ex- 
actly conserve mass, which would not be the case in 
a strictly linearized setting); we have introduced the 
peculiar-volume deformation, following from (38), 

a(t, x''):=i+ am + e mi + e , (42) 

li := I (^'',) ; Hi := II (^",) ; Illi := III (^",) , 

(43) 

where the principal scalar invariants of the matrix .^J^^j 
are given by 

II (^".) := .^%.^\5\ , 

IIl(^'^,) := ieaoce^^'-^",^", 



(44) 



1. The kinematical backreaction 

In addition to the local volume deformation J, we need 
the invariants of the expansion tensor 8 j-. As we have 
introduced a background by Eq. (30), we can define a 
peculiar-expansion tensor 9^^ := — H{t)S''j with re- 
spect to this background. The three principal scalar in- 
variants of the full expansion tensor decompose into in- 
variants of the peculiar-expansion tensor as follows: 

i(e^) = 3i/ + i(0^), 
ii(e*^) = 3H^ + 2H i{e'j) + ii{0'^) , 

111(6'^) ^H^ + HH{e'j) + mi{d'j) + 111(6' . (45) 

Inserting Eqs. (45) into Eq. (11), we can write the back- 
reaction variable in terms of invariants of the peculiar- 
expansion tensor: 

Qv = 2{lW^))^-^{W^))l. (46) 

This nontrivial result demonstrates that the backreaction 
effects do not depend on an assumed homogeneous ref- 
erence background (a Hubble flow): backreaction is only 



due to inhomogeneities. Now, using the formula, Eq. (8), 
the backreaction term (46) can be expressed in terms of 
the formal average (7): 



Qv = 



(47) 



Now we have all we need. Plugging Eq. (38) into Eq. (25) 
and subtracting H{t)S''j, we obtain for the scalar invari- 
ants of the peculiar-expansion tensor. 



RZAj(^z^) 



RZA 



11(9%) 



RZA 



IIK 




6 e(i) 

which implies for the backreaction term 



RZA 



■D 



i 0)2 




mi] 

my ' 

im\ m 
^m) m ' 

(48) 



(49) 



which with Eq. (42) yields a compact form of the RZA 
backreaction model: 



RZA 



Qv 



(71 + ^72 + e'73) 



with: 



71 
72 

73 



= 2(IIi)^-|(lO^ 



6 (Illi) 



(ni>x(ii; 



2(Ii)^(IIIi)^-|(IIi)f . 



(50) 



We note that the initial metric tensor G does now not 
explicitly appear in the expression for the backreaction 
model; the above expression is formally equivalent with 
the Newtonian expression [BKS]. 



2. The intrinsic curvature model 

Unlike in the Newtonian case [BKS], where there 
exists a flat embedding space, RZA is fully intrinsic 
(i.e. the structures are not embedded into an exter- 
nal space, but propagate within the space given by 
the coframe deformation), and we have a non- vanishing 
scalar 3— curvature. Its expression may be found by com- 
bining the Hamilton constraint and the Raychaudhuri 
equation: 



6A , 



(51) 



n = 6H(e^) - 4F(e*j) - 41(9*^; 

or in terms of the first two scalar invariants of the 
peculiar-expansion tensor: 

6fc 



n ^ 611(9]) - 4^(9]) - 4i(0*) - 12Hl{9i) - 



am 



(52) 
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For the averaged peculiar-scalar-curvature, Eq. (13), us- 
ing Eq. (48) and Eq. (52), we get: 



RZA 




0)2 



-6(A-^). (53) 



Inserting the result for 3 from Eq. (42), we find the 
explicit time-dependence in a form similar to Qx> in 
Eq. (50): 



RZA 



(71 + + e'73) 



-6 



with: 



71 

72 := 
73 



-2(IIi)^-12(Ii)^|-4(Ii)^^, 
= -6 (Illi)^ - 24 (Ili)^ f - 8 (11;)^ X , (54) 



-36 (IIIi)^|-12(IIIi)^X 



In this expression there are two important formal dif- 
ferences compared with the functional form of Q-p. First 
of all, the time-derivatives of ^ no longer disappear au- 
tomatically from the solution. In addition, the curvature 
term also explicitly depends on the background via H. 
This leads to a time-dependence of the coefficients 7^, 
whereas for the kinematical backreaction functional they 
were simply constants. 

In the case of the Einstein-de Sitter (EdS) background, 
the growth of the first-order perturbation goes with the 
scale factor q = a, and so we can simplify the expression 
to 



RZA 



Wv = 



As e 



-lOHii^ (ij)^ + 2 (iij)^ e + 3 (iiij) J e 

2(IIi)^ + 6(IIIi)^e 

^ i+^{h)^+e{ih)x+e{mx ^ ' 



i^^ in the EdS case, we recover the well-known 
result that the leading curvature contribution goes as 
cx a~2. xhe second term, which is the leading second 
order contribution, goes as as expected [46]. We will 
learn below that this form of the peculiar-curvature func- 
tional is not the best choice. As it will turn out in the 
explicit consideration within the class of exact spherically 
symmetric solutions. Section VB 1, there exists a better 
approximation for numerical evaluations that employ Qv 
together with the exact integrability constraint. This ex- 
pression and the comparison with the above expression 
will be explicitly provided in Section IV B. 



IV. CONSISTENCY CHECKS 

In order to evaluate the goodness of the averaged model 
based on the RZA deformation, we will perform three 
consistency checks in this section. 

A. Newtonian limit 

First we show that the RZA result of Eq. (50) for the 
backreaction model has the correct Newtonian limit. The 
procedure to arrive at the Newtonian limit of the rela- 
tivistic quantities follows the prescription of [LI], that 
sends the nonintegrable Cartan coframes to integrable 
ones: 



'7" 



N a _ ra 
V i — J I, 



As was demonstrated in [LI], this leads exactly to the 
Newtonian system of fluid equations in the Lagrangian 
description. is the position vector field that maps 
the Lagrangian coordinates to Eulerian positions, and 
its Lagrangian gradient encodes the volume deformation 
of fiuid elements. In ]BKS] this has been approximated 
by the Newtonian form of the Zel'dovich approximation. 
Comparing the corresponding expressions we find that 



N 



^1- 



i.e. the time-derivative of the displacement one-forms 
can now be represented by spatial derivatives of a poten- 
tial -0 in the Newtonian approximation. The invariants 
of Eq. (43) similarly become 



li :=I(/;^.), Hi :=II(0'|;.), Illi :=iii(^'j;.) 



(56) 



With this spatial geometrical limit we also send the gen- 
eral curved space, on which we defined our average, to 
Euclidian space. Note that the existence of a vector field 
/" implies that the counter-index a becomes a coordinate 
index of (now existing) global coordinates in an embed- 
ding space. The Riemannian volume average automati- 
cally corresponds to the Euclidian volume average over a 
flat domain Vi. (Note that sending the speed of light to 
inflnity is only needed in order to change the spacetime 
metric signature; since the backreaction model is spatial, 
and since we have eliminated all terms that would con- 
tain the speed of light (the curvature), we have not to 
send c to infinity.) 

Taking the spatial geometrical limit is thus sufficient to 
reduce the result, Eq. (50), to the Newtonian expression: 



NZA 



with: 



(l+C(Ii>i,,+e (Hi) 



e (Illi) 



Ti := 2 (Ili)^^ - I (li)^^ , 

T2 :== 6 (Illi)^^ - f (li)^^ (Ili)^, , 



(57) 



T3 := 2 (li)^^ (Illi)^, - I (Hi 
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B. Matching curvature expressions 

Another consistency check considers the scalar curva- 
ture in the RZA. There are in principle three ways to 
calculate the intrinsic curvature that all agree for an ex- 
act solution. The first way is to use the RZA metric and 
to calculate the curvature geometrically. The other two 
ways will be studied below; both relate the scalar cur- 
vature to extrinsic curvature invariants (note that these 
are relevant for this work, since we study kinematical 
backreaction being related to extrinsic curvature): 

a) using a combination of the Hamilton constraint and 
the Raychaudhuri equation, 

b) using the integrability condition. 

Thus the scalar curvature is a good tool to test the RZA: 
we expect that to the first order the two results obtained 
by a) and b) are the same. The expression for the way 
a) has already been obtained as Eq. (53). For the way 
&) we express Wv through Eq. (14) as: 



(58) 

Using (33), (49), the Friedmann equations for the back- 
ground scale factor, and 



RZA 



ait) 



(59) 



one can show after a long but not too technical calcula- 
tion the following: 



RZA 



Off 



+3 U4^ 




-6 1 ^ 

a. 



1 



dt' 



(60) 



Using this expression for the averaged curvature assures 
the integrability condition (12) to hold for a given kine- 
matical backreaction model. In other words, this expres- 
sion respects the conservation law for the combined ac- 
tion of kinematical backreaction and averaged curvature 
(third equation of Eq. (17); note that the averaged rest- 
mass is individually conserved). 

If we now compare Eq. (53) and Eq. (60), we find that 
the RZA is consistent if and only if the remaining inte- 
gral vanishes. After a change of integration variables this 
integral has the form: 



6 a{i) (y)f 



d2 



. (61) 



As it is zero for all values of already the integrand 
has to vanish. Its first three factors are nonzero which 



(62) 



Using the definition of (J)^- this implies that the RZA 
approximation is consistent iff: 




= 1 



= (Illi)^ 



(63) 



which encodes the fact that it is strictly speaking only 
a first-order scheme. Its success in the Newtonian case, 
however, motivates its use in the relativistic case as well, 
despite this result which otherwise stated reads: 



q{u) (iii)x + 3e(nii) 



«(0 



0)f 



'-d^ = O m)^) . (64) 



This shows that the deviation from Eq. (53) contains only 
terms of second and higher order in the initial conditions. 

An equivalent way to check the consistency is to 
insert ^^^Qv (Eq. (50)), ^^^Wv (Eq. (54)) and 
^^^^©(Eq. (59)) directly into the integrability condition 
Eq. (12). The result may be written as 



6go»„(l + goO ■ d^ 



0), 



. 



This can be interpreted as the amount by which the "clo- 
sure" of the integrability condition of Q-p and ^ "^^Wv 
fails. As the prefactors are nonvanishing, we recover the 
condition of Eq. (62). 



C. Self— consistency test in terms of the scale factor 

We have already seen above that the RZA delivers 
only consistent expressions, if it is employed strictly as a 
first-order scheme. Therefore, an extrapolation as sug- 
gested in [LI], bearing a number of advantages as dis- 
cussed there, is only an approximation for the terms at 
second and higher order. It is therefore not surprising 
that, in addition to the different curvature expressions, 
we also have two concurring definitions for the volume 
scale factor. One is the kinematical volume scale factor 
defined as the cubic root of the local volume deformation 



in the RZA, Eq. (41), ^^^av = (^^^^J)^^^ The other 
one is ^^^au, calculated from the RZA of the backreac- 
tion term and using this latter as a source of the gen- 
eral equations governing the average evolution, i.e. the 
prescription of Subsect. VIA. Just like the curvature ex- 
pressions, they do not coincide for the RZA. For a full 
n"^-order calculation they are the same up to the given 
order. For the curvature we will see in Subsect. V B 1 
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that the expression derived from the backreaction term 
is more powerful and seems more appropriate to capture 
the nonhnear evolution. With this insight we are going 
to use ^^^ap in the concrete calculations presented fur- 
ther below. Fortunately, it will turn out that the possible 
error induced by instead calculating the kinematical vol- 
ume scale factor is numerically rather small: evaluating 
the quantity 



IRZA 



KIN 



RZA 



(65) 



we found that ea stayed well below 0.1 at all times, as 
long as ^'^^aj} did not approach zero and the domain's 
effective radius was larger than 16 Mpc. 



V. EXACT SUBCASES 

We are now going to study subclasses of exact solutions 
that are contained in the above approximation. We will 
learn that (i) a subclass of the averaged plane-symmetric 
solutions, (ii) a subclass of the averaged LTB solutions, as 
well as (iii) a subclass of the averaged Szekeres solutions 
are found by suitably restricting the initial data. This 
fact demonstrates that the investigated approximation 
interpolates between two kinematically orthogonal exact 
solutions that are, however, unlike in the corresponding 
Newtonian model, subjected to constraints. This latter 
fact suggests that one may be able to construct more 
general backreaction models that include these cases in 
full generality. However, we do not expect this to be 
possible for first-order deformations. 



A. The 'Relativistic Zel'dovich Approximation' 
and plane collapse models 

In Newtonian cosmology the 'Zel'dovich approxima- 
tion' is an exact three-dimensional solution to the Newto- 
nian dynamics of self-gravitating dust-matter for initial 
conditions with ll{9ij) = = 111(0^^) at each trajec- 
tory [9] . This " locally one-dimensional" class of motions 
contains as a subcase the globally plane-symmetric so- 
lution. In the relativistic case one may study the cor- 
responding model using the plane-symmetric ansatz for 
the line-element, 



-dt^ + a {tf (dx^ + dj/2 + (1 + p (z, t)f dz^^ 



(66) 

which also has vanishing higher invariants of the pecu- 
liar expansion tensor: II (6*^) = = lll{dij). The first 
invariant is nontrivial and reads 



I e 



l + P{z,t) 



(67) 



The equation determining the time-evolution of P [z,t) 
was in the Newtonian case simply 

e + ef = -a^Gq + a , (68) 



which gave 



P (z, t) + 2-P (z, t) = AnGQH {P {z, t) - Po (z)) . (69) 

Hence, the Newtonian plane collapse had two solutions, 
e.g. for an EdS background, 

P{z,t)^Po{z) + aCi{z) + ^^ , (70) 

a growing and a decaying one. In the relativistic case, 
however, there are more constraints. In the Lagrange- 
Einstein-System of [LI], also a link to the scalar curva- 
ture comes in: 



6' + ee; ^ (AirGg + A) (5; - K) 



- Qfei = 167rGg + 2A - 7^ 



(71) 



(72) 



These two equations combined also give Eq. (68) for the 
relativistic case. Additionally, however, they have to be 
satisfied individually. As the plane-symmetric metric 
ansatz Eq. (66) implies that TZ'^j = 0, the relativistic so- 
lution space is not the same as the Newtonian one. For 
the Hamilton constraint Eq. (72) we find 

-P (z, t) = -AirGgn {P (z, t) - Pq (z)) , (73) 
a 

which is now, for 7?. = 0, only a differential equation of 
first order with the solution 



P{z,t)^Po{z) + 



Cjz) 

.3/2 



(74) 



This P{z,t) also satisfies Eq. (71), but one of the solu- 
tions of the Newtonian case has disappeared. 

From (74) it follows that P{z,t) = Po (z) + q{t)C{z), 
with q being a solution of (33). Inserting this result into 
(67), and using Eq. (46), we find for the backreaction 
term 



gplanc^_ - 



2/ 



V 3 v 1 + e (ii>x 



(75) 



where we recover ^ (t) ~ {q (<) — q (to)) /q (to)- Siiid where 
Ii = l{9 j){ti). Extracting the coframe from the met- 
ric (66) we find also that Ii = I(^°(z,<i)j, where 



{z,ti) is the RZA perturbation. This shows that the 
plane-symmetric metric is, as in the Newtonian case, a 
particular exact solution that is contained in the solu- 
tions of the RZA. Note again, however, that this solution 
in the RZA as well as for the plane-symmetric metric 
does not have the growing mode that was present in the 
Newtonian solution. This is due to the vanishing scalar 
curvature for cylindrical symmetry and the relation to 
the Hamilton constraint that did not exist in the New- 
tonian case. (Note that the integrability condition (12) 
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is satisfied by this solution.) This is another interesting 
example of a case in which a class of Newtonian solutions 
may not automatically provide a solution of general rel- 
ativity. 

For negative (li)^;, corresponding to over-dense re- 
gions (Eq. (92)). 2^'*^"° is diverging at some time when 
1 + (,{t) approaches zero, even though the solu- 
tion is decaying. Our special initial conditions imply 
a one-dimensional symmetry of inhomogeneities on a 
three-dimensional background (cylindrical symmetry), 
and the diverging Q^'''"'^ is supposed to mimic the highly 
anisotropic pancake collapse in the three-dimensional sit- 
uation. 

The plane-symmetric case has also been considered in 
[63] , but in the framework of a post-Newtonian approx- 
imation. The evolution model discussed here does not 
need to invoke a post-Newtonian approximation, but this 
can be useful for a precise calculation of the initial invari- 
ants. 



The 'Relativistic Zel'dovich Approximation' and 
spherical collapse models 



the averaged scalar invariants of the expansion tensor can 
be calculated: 



Ait 



dr RR^ 



J'/LTB Vltb Jo VTT2E 
dr ( R'^R 

(II(0V)ltb 



An 

VhTB Jo 



(iii(e^)) 



47r 



VTT2E 
. dr(R^ 



dr ; (80) 
dr ; (81) 
:dr , (82) 



J'/LTB s^LTB Jo VT+2E 
where the Riemannian volume of an LTB-domain is given 

by 



3 Jo y/l + 2E 



(83) 



The deviation from constant-curvature can also be aver- 
aged on a LTB-domain: 



Wltb = 



LTB 



2/3 



(84) 



LTB 



1. The LTB solutions 

We now investigate the Lemaitre-Tolman-Bondi 
(henceforth LTB) solutions [45]. The LTB model, written 
below in the time-synchronous metric form, is the gen- 
eral inhomogeneous spherically symmetric solution for ir- 
rotational dust. The spherical solution can be seen as a 
superposition of infinitesimally thin homogeneous shells 
governed by their own dynamics. In such a domain one 
can show (see ]29] for a demonstration but with different 
notation, and the review [62[ for a comprehensive study of 
backreaction in the LTB solutions) that the line-element 
has the form: 

= -d^' + rSlrTdr^ + ^'(^'0d^^^ (76) 

1 -I- 2h[r) 

E being a free function of r satisfying E{r) > —1/2; the 
prime denotes partial differentiation with respect to r. 
In this metric, the scalar parts of Einstein's field equa- 
tions read: 



with 



ATrp(t, r) 



M'{r) 



R'{t,r)R'^{t,r)' 



and 



1-2, , GM(r) 



(77) 



(78) 



M being another free function of r; the dot denotes par- 
tial time-derivative. Using the relation between the co- 
efficients of the expansion tensor and the metric tensor. 



6' 



1 



(79) 



(7^) 



LTB 



167r 
Vltb 



dr [ER) 



dr 



(85) 



There are two cases for which we find relations between 
the invariants without having to solve the system (77) 
and (78) explicitly: the first one is a separable solution 
R {t, r) of the form 

R{t,r)^f{t)-g{r) . 

The second one is the case of an LTB domain where 
E{r) = E. The restriction E = const, corresponds to 
self-similar LTB solutions if we require at the same time 
that the function A/(r) oc r [61[. 

In both cases, one can show for R {t, 0) ~ 0: 



(^m)kTB - 1 (i(0^.))Lb 

("I(0'.))ltb = ^(I(0'.))ltb 



(86) 



Combining these terms in the backreaction Qltb given 
by Eq. (11), we get for a spherically symmetric E = 
const. -domain or a separable R (t, r): 



Qltb = , Wltb = 



(87) 



where the result Wltb = follows from Qltb = by 
the integrability condition (14) and its definition (13). 
We here generalize a result obtained in [51[ to a non-flat 
domain for a special case. 

Note that, using the expression for the curvature of 
Eq. (54), we find that the result is not zero for the invari- 
ants Eq. (86). This means that this expression is not yet 
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containing the correct second-order contributions. This 
is one important reason for why we use in the evalua- 
tion of the importance of Wv in Section VI B 2 the nu- 
merically integrated expression, starting from the Q-p of 
Eq. (50). 

In the case E — const., Eq. (84) and (13) give for the 
averaged scalar curvature 



(7^) 



LTB 



12£: 



(88) 



With Wltb = 0, one can express fcp. as a function of E: 
2E 1 



k-D, = - 



< 



(89) 



As i? is a growing function of rx), kxi^ becomes smaller 
when we increase the averaging domain. 



2. The Szekeres solutions 

A possible generalization of the LTB models is 
the quasi-spherical Szekeres model that has additional 
anisotropics which destroy the spherical symmetry of the 
LTB models [5]. The line-element reads: 

ds^ = -dt^ + X^dr^ + {dx^ + dy^) , (90) 

and with = x + ty ; Q = x — ly, we have 

R{t,r) 



£(r,C,C) =aWCC + 6(r)C + c(r)C + rf(r) ; e = 0,±l, 

where a prime denotes dr- We consider only the quasi- 
spherical model e = 1 in which case, for £' = 0, it simply 
reproduces the LTB model. But even for the general 
case £' 7^ it has been shown in [5] that the invariants 
take the form of Eqs. (80) to (82) with k (r) = -2E (r). 
Therefore also the flat quasi-spherical Szekeres model has 
no backreaction; 



Qsz = , Wsz = 



(91) 



We considered only this example for the flat case, but the 
other flat cases should be accordingly solvable. However, 
in the really interesting cases, i.e. for arbitrary curvature, 
it is highly difficult to find a general statement, as it 
is also the case for the LTB model. For further work 
on the Szekeres model see [42, 48], and in relation to 
observations [34, 35, 50[. 



the averaged equations it then follows that, e.g. a flat 
spherically symmetric but inhomogeneous domain, cut 
out of a homogeneous FLRW model, has no influence 
on the kinematical expansion properties. Such a result 
is known as Newton's 'Iron Sphere Theorem' in a New- 
tonian model. Let us now look at situations where the 
backreaction term of the RZA approximation vanishes. 
Using it in the form of Eq. (50), one can show the 
following propositions: 



Proposition 1. 



RZA 



(ni)x = ^(ii)i 



(nii) 



Proof, see Appendix A. 



Proposition 2. 



RZA 



(""^n(e',))^ = ir^i(ev)); 



(H-^iii(e',))^ = ^rAi(e',)); 



Proof, see Appendix A. 

This means that, for a flat background, the RZA cor- 
rectly reproduces the average evolution of the LTB mat- 
ter shells, without knowing the specific distribution of 
matter M (r). This corresponds to the 'Iron Sphere The- 
orem' of Newton in the relativistic case (compare also 
[51[), and it is in the spirit of Birkhoff's and almost- 
Birkhoff theorems, but for non-vacuum spacetimes (see 
the recent papers [32, 33[ and references therein). 

In [BKS[ it was shown that, in the Newtonian case, 
the averaged Zel'dovich approximation also describes the 
behavior of a spherical matter distribution. Here we find 
the analogous property for the flat LTB and flat Szekeres 
metrics. 

To summarize one may say that the RZA approxima- 
tion interpolates between the exact GR solutions of a 
plane-symmetric metric, the flat LTB and the flat Szek- 
eres metrics, which adds reliability to the employed ex- 
trapolation of a strictly first-order scheme and its use 
within the exact averaging framework. 



VI. THE EVOLUTION OF COSMOLOGICAL 
PARAMETERS 



3. GR theorems corresponding to Newton's 'Iron Sphere 
Theorem ' 

An interesting property arises, if we look at situations 
where the kinematical backreaction term vanishes. From 



The backreaction term itself decays in the averaged 
"Relativistic Zel'dovich Approximation". However, to 
quantify the deviations from the behavior of the scale 
factor and especially its time-derivatives in the standard 
model, the different strengths of the sources in the gen- 
eralized Friedmann equation have to be compared. As 
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the matter source term decays faster than the backreac- 
tion term, the influence of this latter grows. To evalu- 
ate the importance of this growth in the standard cos- 
niological picture, starting with a nearly homogeneous 
and isotropic initial state, we only need to determine 
the magnitude of the three invariants of the perturba- 
tion one-forms The time-evolution of ^^"^Qd is 
then determined by Eq. (50). It then also determines the 
evolution of all other cosmic parameters via the averages 
of the Hamilton constraint, Eq. (10), and Raychaudhuri's 
equation, Eq. (9). The context of the standard scenario 
implies that the calculation of the initial values is per- 
formed in a universe model that is close to spatially flat, 
looking back to a history that is identical to the stan- 
dard model. If we additionally neglect tensor modes, this 
means that we are in the limit described in Sect. IV A. 
We therefore assume that we can use the values of (Ij) p. , 
(Ili)p. and (Illi)j,. in [BKS] as our initial conditions for 

the averaged invariants of to a very good approxima- 
tion. In formal analogy of Eq. (50) to the corresponding 
expression of [BKS], this then implies that many results 
of [BKS] carry over to the RZA context. There are, how- 
ever, new phenomena that emerge due to the fact that, 
unlike in the Newtonian ZA of [BKS[, geometry is a dy- 
namical variable and the RZA develops non-vanishing 
scalar curvature. In this section, we will therefore com- 
ment on which results of [BKS[ remain valid and discuss 
where the GR description brings in new phenomena. As 
we emphasized already in the introduction, an emerging 
curvature is key to the new interpretation of the dark 
components in the standard model. 



A. Calculation of the scale factor a-o 

We calculate the average scale factor ^^^au as in 
[BKS] directly by integrating the averaged Raychaudhuri 
Equation Eq. (9). The input is the average density 



'\Q)v = -rSH (1 + {S (ii))i,J = -^QH (1 - iWv,) 



V 



(92) 

and the RZA backreaction model, Eq. (50). The ini- 
tial conditions for ^^^a-D are ap- = 1 and dp {ti) = 
d(ti)(l + i(Ii)^J. 



positive or negative. These invariants fluctuate with the 
variance, e.g., cr^{R) = E (Ii)x)^ • In our calculation of 
the time-evolution of axi{t), we consider one-cr fluctua- 
tions of the averaged invariants for spherical domains of 
radius R, e.g. (li)D^ = ±cri (R). ]BKS] showed how these 
fluctuations are linked to the matter power spectrum. As 
indicated, (7i{R) will explicitly depend on the radius of 
the initial domain, but implicitly also on the shape of 
the power spectrum. [BKS[ used a standard CDM power 
spectrum normalized to (Tg = 1 and an h of h ~ 0.5. 
In addition to the EdS background considered there, we 
will also present some results for a standard ACDM back- 
ground with ~ 0.73, h = 0.7 and CTg = 0.8. 

Although the possible initial conditions have a rich 
structure due to the fact that we can choose any com- 
bination of the signs of the three initial invariants, we 
choose for the presentation equal signs for all the invari- 
ants, i.e. for expanding, under-dense domains, ai > 0, 
all initial invariants are taken to be positive, and for col- 
lapsing, over-dense domains, ai < 0, all initial invariants 
are taken to be negative. (A show-case of mixed signs 
has been analyzed in [BKS[.) 

It is important to recognize that our choice of using 
one-cr fluctuations in the invariants does not imply that 
the fluctuations of other parameters constructed from 
them are then also one-cr fluctuations of these param- 
eters. For the initial value of Q-p for example a one-cr 
fluctuation of the invariants is related to a one-cr fluctu- 
ation of Qd by 



(ii>;, (Hi) 



8 



(94) 

On the other hand, calculating Qp. with one-cr fluc- 
tuations yields 



RZA 



Qv, = 2au {R) - -<^! (R) . 



(95) 



In an abuse of language we will nevertheless speak of one- 
cr fluctuations in the following, but one should keep in 
mind that we mean those of the initial invariants (except 
for Fig. 1). For the values of ^^^Qvi: in the case of large 
values of i?, the two prescriptions coincide as crj (i?) drops 
off faster than au (i?) . 



1. Statistics of initial conditions 

As discussed, we use for the numerical evaluation the 
initial values given in [BKS[. As was shown there, the 
expected values of the initial invariants are zero, 

E [(li)^J = E [(Ili)^J = E [(Illi)^J = 0. (93) 

E [. . . ] denotes the ensemble expectation value over many 
realizations of universe models. However, for a speciflc 
domain, any of the volume-averaged invariants may be 



2. Scale-dependence of initial conditions 

From the very definition of the averaged parameters 
it seems natural that all averaged quantities would be 
scale-dependent and decay with a growing domain V. 
That this is not necessarily the case shows the example 
of the initial expansion- and shear-fluctuations. In the 
RZA they are given by 



14 



Using the acceleration equation in terms of coframes of 

[LI], one can show that for the RZA we have 1 = 

—(5, where 6 is the local density contrast. This means that 
also in the RZA, as already in the Newtonian case, the 
expectation values of expansion- and shear-fluctuations 
are no longer scale-dependent. They are rather given by 

(e^)^ {0)1] = (^^^ d'kn (k) - al (i?)) ; 

(98) 

E = I d'fc^i {k) , (99) 

where we again made use of the assumed approximate 
flatness of the initial universe model, necessary to em- 
ploy the Fourier transformation. Interestingly, only part 
of the expected expansion fluctuations still contains the 
information on the domain T). The other part and the 
shear is domain-independent. Calculating the value of 
this integral may be used to estimate the importance of 
backreaction, if the shear fluctuations were negligible. To 
this end we calculate 



E 



Q.trunc 



1 



9H^ 



(100) 



where stands for today and we evolved the model from 
the initial time up to today with the leading mode 
of Eq. (50) only. A quantitative estimate with an expo- 
nential IR cut off at the Hubble scale and UV cut off 



at Ikpc then yields 



0.73. This il- 



Q.trunc 

lustrates the well-known fact that the expansion- and 
shear-fluctuations by themselves are important even in a 
perturbative framework. In the backreaction term, how- 
ever, they combine in a way that leaves only the domain- 
dependent contribution in the second term of the expan- 
sion fluctuations. In the Newtonian framework this is 
expected by the fact that Qd can be written as a surface 
term. In GR it is not a necessity, but in the linear RZA 
the cancellation is still effective. For higher orders, how- 
ever, this is no longer true and [23] reported the survival 
of a domain-independent contribution to at second 
order. 

To close this section we calculate, as an illustration of 
the scale-dependence of the parameters, the backreac- 
tion term Qp. It will imprint its scale-dependence on 
the other parameters and is therefore particularly inter- 
esting. To get a feeling for the magnitude of the values, 
we evolve it again with the mode until today and 
normalize it with to get ^Iq' ■ The result is shown 
in Fig. 1. We plot the one-cr fluctuation of fig" with the 
correct cr— interpretation of Eq. (50). The result shows 
that only below the assumed "almost homogeneity scale" 
of about 400 Mpc the (with a~-^ upscaled) backreaction 
term begins to enter the range of a percent contribution. 



■S 4 



EdS, h = 0.5 
ACDM, h = 0.7 



100 200 300 

scale 2R in Mpc 



400 



Figure 1: One-a fluctuations of f^g" as a function of scale for 
(i) an EdS background {h = 0.5, erg = 1) and (ii) for a ACDM 
background (^a ^ 0.73, h = 0.7, = 0.8). 



For larger scales it is clearly below this value, while for 
smaller scales the nonlinear terms in the backreaction 
functional count in. This explains that in the following 
all parameters converge to their background values for 
large I). It is here where our model faces its strongest 
restriction, since we neglect any interaction of structure 
with the assumed background model. 



B. Time— evolution of cosmic parameters 

Having discussed in the previous section, how we fix 
our initial conditions and determine the scale factor from 
them, we will now present and discuss the results. 



1. Evolution of the scale factor a-p, the Hubble- and 
deceleration-parameters 

As mentioned, the evolution of the dynamical quanti- 
ties a-p, Hx> = hxilax) and qx> = —{dv/av)/H^ turns out 
to coincide with those of [BKS[ for our choice of initial 
conditions. This means that the partially drastic devi- 
ation from the background values for these quantities, 
that has been found in [BKS[, also occurs in the RZA. 
In the case of the volume deceleration parameter g-p for 
one-cr fluctuations this leads to deviations of 30% on a 
scale of 200 Mpc. 

The effects on and are smaller, but may also 
become important for special regions, that are more than 
one-cr away from the background. In the relativistic 
framework the deviations of the volume scale factor av 
from the background scale factor a{t) can also be inter- 
preted as giving the strength of metrical perturbations, 
since the calculation of the volume just involves the met- 
ric determinant, and not higher derivatives of the met- 
ric. Taking the cube of this deviation gives us a typical 
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Figure 2: Evolution of tlie volume scale factor, normalized by the background scale factor, on typical domains of 50 and 200 
Mpc, and on a domain that is a 2— sigma fluctuation in the initial conditions on 200 Mpc for comparison. Left: EdS background 
with = 1 (/i = 0.5, as = !)■ Right: ACDM background with f7„i = 0.27 {h = 0.7, ag = 0.8). 
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Figure 3: Evolution of the domain-dependent cosmological parameters of Eq. (18) with cosmic time. One background is the 
EdS model with Qm ^ I (h = 0.5, as = 1) (left), the other one the ACDM model with Qm = 0.27 (h = 0.7, as = 0.8) (Right: 
the background density parameter is plotted here as the upper curve). The figure shows values for an expanding under-dense 
domain of 200 Mpc effective diameter with one-o" fluctuations of the initial invariants of the perturbation one-form. This figure 
confirms the findings of [BKS] reporting substantial deviations from the background of e.g. the matter density parameter, while 
the quantitative importance of the backreaction parameter is seemingly negligible. The new interpretation in the GR context 
is mirrored by the curvature density parameter that has to compensate (in the case of an EdS background without a 'Dark 
Energy' component) the large deviations in the matter density parameter. 



strength for the volume fluctuation, e.g. for one-sigma 
fluctuations on the scale of 200 Mpc we find a 12 percent 
effect, see Fig. 2 for the evolution of the scale factors. 
This also means that the influence of metrical perturba- 
tions is not necessarily small, if their averaged effect is 
considered (see the related discussions in [21, 27, 55]). 
However, it is still subdominant compared with the over- 
all backreaction effect that depends on second derivatives 
of the metric. 

In the next subsection we will discuss the results for 
the density parameters that are relevant for the inter- 
pretation of the energy budget of the Universe and its 
relation to the 'dark components' in the standard model. 



Recall that we quantify fluctuations on Lagrangian do- 
mains in this paper; for a quantification of fiuctuations of 
cosmic parameters on domains that correspond to actual 
observational geometries see [65]. 

2. Evolution of the density parameters 

Let us consider the cosmological density parameters, 
Eq. (18), describing the energy content balance of the 
universe model. In the average scenario they are domain- 
dependent quantities. In Fig. 3 we show their evolution 
with cosmic time on an expanding typical (i.e. one- 
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Figure 4: This figure corresponds to Fig. 3, but shows the corresponding values for an expanding domain of two-cr fluctuations 
of the initial invariants of the perturbation one-form. 
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Figure 5: Evolution of the domain-dependent cosmological parameters of Eq. (18) with cosmic time. One background is the 
EdS model with fi™ = 1 (/i = 0.5, ag = 1) (left), the other one the ACDM model with fi^ = 0.27 [h = 0.7, erg = 0.8) (Right: 
the background density parameter is plotted here as the lower curve). The figure shows values for a collapsing over-dense 
domain of 100 Mpc effective diameter with one-o" fluctuations of the initial invariants of the perturbation one-form. 
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Figure 6: This figure corresponds to Fig. 5, but shows the corresponding values for a collapsing domain of one-cr fluctuations 
on the scale of 50 Mpc. On this scale we appreciate a singular pancake collapse for the EdS model; Fig. 7 illustrates that the 
backreaction term now becomes not only qualitatively but also quantitatively significant. 
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Figure 7: Ratio of tlie cosmic parameters for backreaction and curvature to the matter density parameter for a scale of 50 Mpc 
as plotted in Fig. 6. As during the pancake collapse Ht> becomes 0, the parameters themselves are no longer well-defined. In 
the ratios plotted in this figure, however, cancels. 



sigma) domain 200 Mpc. Fig. 4 shows the situation again 
for an expanding domain of 200 Mpc, but here found with 
two-sigma probability. In Figs. 5 and 6 we give examples 
of smaller typical domains of 100 and 50 Mpc, but this 
time in the collapsing phase. 

For the case of the EdS background with Ha = = 
and rjfc = 0, the global matter density parameter is Vim = 
1. Regionally, however, the plot shows that fi^ is on 
200 Mpc domains typically lower by 20%, i.e. this scale 
is dominated by under-dense voids. In the relativistic 
framework the lower matter density is compensated by 
an emerging curvature parameter f2^. 

It is interesting that the curvature deviation from the 
background curvature, Wu, a second-order quantity, can 
be a lot bigger than the second-order quantity Qt>. This 
is due to the fact that in the perturbative expansion the 
second-order terms both scale as a~^, and in view of the 
integrability condition Eq. (12), the second-order contri- 
bution to Wv is a factor of 5 larger than the second-order 
term of Q-p. Therefore, even a small backreaction contri- 
bution of only 2% will already lead to a 10% modification 
of the averaged curvature. 

For the ACDM background the matter density param- 
eter is also reduced with respect to the background value 
and again we have curvature emerging from a flat back- 
ground. Today, however, it is not as big as in the EdS 
case, since the cosmological constant dominates. 

In both cases we see that, even though the backreaction 
contribution stays tiny (as was already found in [BKS]), 
we have a considerable amount of curvature. Compar- 
ing this curvature contribution with the flat geometry 
of the background on a given scale allows its interpreta- 
tion in terms of 'Dark Energy': the standard interpre- 
tation is that the matter distribution evolves on a flat 
geometry; hence one would add a fundamental compo- 
nent to compensate the actually existing curvature that 
we model here. In other words, the matter distribution 
has to be seen on a curved space section and not on a 



flat background. Taking this point of view, the emerg- 
ing curvature quantifies, on a given scale, the amount of 
'Dark Energy' that would be needed to compensate it in 
a quasi-Newtonian model. We find for the example of 
Fig. 3 that our model predicts the existence of typical 
domains with a diameter of 200 Mpc of about 1 /4 of the 
needed amount of 'Dark Energy'. This amount can in- 
crease to 0, 4, if the 200 Mpc domain is slightly untypical 
(i.e. found with a two-sigma probability. Fig. 4). 

Figs. 5; 6 and 7 show that the backreaction can have 
the opposite effect by looking at smaller scales. A typ- 
ical, collapsing domain produces on average a positive 
curvature (corresponding to a negative curvature density 
parameter) . In light of the interpretation above, the col- 
lapse produces a large amount of 'Dark Matter' in the 
form of positive curvature. While the 'cosmological pa- 
rameters' start to loose their sense on this scale, we can 
clearly see, e.g. in Fig. 5, that the curvature contribu- 
tion is of the order of the density contribution, since it 
compensates the produced over-densities. The physical 
interpretation is also clear: over-densities are Irosted in 
positive-curvature environments. The order of magni- 
tude of 'Dark Matter' is also comparable with that of 
'Dark Energy' in expanding domains, i.e. we can also 
here say that a substantial fraction of the density pa- 
rameters is contained in the curvature parameter, while 
it is of comparable magnitude and not dominating. 

Given the conservative assumptions of our model these 
values point to a highly significant effect of backreaction. 
Except in particular LTB models, where 'Dark Energy' 
can be fully replaced, the only generic model that car- 
ries such a large effect has been investigated by Enqvist 
et al. [30] by employing gradient expansion techniques 
that allow to go substantially beyond standard (Eule- 
rian) perturbation methods. These latter have to assume 
a universe model that stays close to the assumed back- 
ground. For comparison, we wish to point the reader 
to a figure in [30]: Fig. 8 corresponds to their Fig. 4 and 
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Figure 8: Scale-dependence of the X— Matter component, := Qq + Q,^, today for typical (one-sigma) domains in the 
averaged RZA. Left: under-dense regions. Right: over-dense regions. As in this case the Hubble parameter goes through zero 
for small scales, we plot the ratio of the X— Matter to the matter density parameter. In this ratio H-d drops out. 



shows the scale-dependence of the X— Matter component 
n'^ := 17^ + today: the X-Matter produces 'Dark 
Energy' on expanding, under-dense domains (left panel) , 
starts to compete with the matter-density on scales be- 
low 100 Mpc and it is still significant (6%) on domains of 
400 Mpc; it produces 'Dark Matter' on collapsing, over- 
dense domains (right panel), starts to compete with the 
matter-density again below the scales of 100 Mpc, and 
shows similar significance as 'Dark Energy' on domains 
with diameter 400 Mpc. Note that the large-scale value 
is astonishingly big, since we are looking at a scale that 
is considered to comply with homogeneity. 



VII. RESULTS, CONCLUSIONS AND 
OUTLOOK 

Based on a nonlinear extrapolation of a first-order rel- 
ativistic Lagrangian perturbation scheme investigated in 
[LI], we modeled the fiuctuations in extrinsic curvature 
in the form of the kinematical backreaction term Qx> for 
the description of the average properties of irrotational 
dust models. We provided backreaction and intrinsic cur- 
vature expressions as nonlinear functionals of the first- 
order Lagrangian deformation and used them as an input 
for the general framework for the average dynamics. The 
large-scale behavior of the backreaction variables, Q-p 
and Wp, at second order is identical to the leading mode 
in the second-order perturbation theory (proportional to 
for an EdS background (while the leading curvature 
contribution is proportional to a^- in conformity with the 
perturbative calculations of [46]). We showed that the 
backreaction model contains in limiting cases the New- 
tonian approximation investigated in ]BKS], as well as 
special classes of exact GR solutions. We argued that 
this backreaction model is a powerful approximation due 
to the successes of the corresponding elements of this ap- 
proximation in Newtonian cosmology (like the large-scale 



performance of the Newtonian Zel'dovich approximation 
in comparison with N-body simulations), and also due to 
the above-mentioned property of interpolating between 
exact GR solutions with orthogonal symmetries. 

We discussed how substantial results of ]BKS] can be 
translated to the relativistic context, the new feature be- 
ing an emerging intrinsic averaged curvature. This trans- 
lation is possible due to the fact that the Lagrangian de- 
scription offered in ]L1] features a clear-cut Newtonian 
limit by sending the coefficients of the coframes to their 
Euclidian counterparts, 77" — !• /°^, where the existence 
of the vector field /° in the Newtonian case describes 
the embedding into a global Euclidian vector space. The 
so-called electric part of the Lagrange- Einstein system 
]L1] is formally identical to the Newtonian equations 
in Lagrangian form, where in the GR case the non- 
integrability of the coframes is responsible for the emerg- 
ing curvature. 

We quantified the fiuctuations of various cosmologi- 
cal parameters on different spatial domains for standard 
model backgrounds (CDM and ACDM power spectrum), 
and we now focus on these quantitative results. 



A. Results 

In order to discuss the quantitative results of our in- 
vestigation, let us concentrate on the averaged scalar cur- 
vature {TZ)t>- The principal idea is to compare the en- 
ergy balance conditions, encoded in cosmological param- 
eters (i) in the standard model and (ii) in the averaged 
model. The balance condition for (i) is entirely defined 
through the chosen background cosmology and furnished 
by the normalized Hamilton constraint. For the two cho- 
sen backgrounds we simply have 

n^ + iiA = 1 , (101) 
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where initially Clm dominates and remains equal to 1 for 
the EdS background, while the importance of A increases 
in the late stages in the background of the 'concordance 
model'. This balance has to be compared with the (scale- 
dependent) balance of the averaged cosmology furnished 
by the normalized averaged Hamilton constraint, 

+ 17? = 1 , (102) 

where il? = + comprises the emerging backre- 
action and curvature due to structure formation. This 
comparison allows us to interpret the X— Matter energy 
density content as a candidate, on some scale, for the 
contribution by that is added ad hoc in the standard 
model. One subtlety of this interpretation is related to 
the scale-dependence of the averaged balance condition 
and also to its evaluation on a chosen background cos- 
mology that is already supposed to contain a substantial 
amount of 'Dark Matter' (for both backgrounds) and, ad- 
ditionally, 'Dark Energy' (for the 'concordance model'). 
We will now discuss the values obtained at the present 
day, and denote the domain by Vq. 

We found, in accord with previous analyses that were 
e.g. summarized in the review papers [16, 21, 24, 27, 
28, 41, 54], that on large under-dense domains of the 
order of 100 — 400 Mpc, fluctuations feature a negative 
averaged curvature due to the fact that on these scales 
the universe model is void-dominated. Note that the 
curvature does not individually obey a conservation law 
like the density; only a combination of averaged curvature 
and backreaction is conserved [17]. 

Let us look for example at a scale of 200 Mpc. While 
the contribution by kinematical backreaction Jig" re- 
mains small in this situation, the curvature density pa- 
rameter (that is positive for negative curvature), fi^" is, 
for one-sigma initial fluctuations, of the order of (e.g. on 
the EdS background) 20% and compensates the lowered 
matter density parameter on this scale, fi^" = 80%. This 
result, being physically plausible, nevertheless paints a 
completely different picture as the one advocated in stan- 
dard cosmology, e.g. [36], claiming that intrinsic cur- 
vature is irrelevant down to the scales of neutron stars 
(see discussions and estimates in ]19, 21, 25, 26, 39, 53]). 
Since this latter prejudice led to the inclusion of a large 
amount of 'Dark Energy' into the background, as mod- 
eled by a cosmological constant in the ACDM model, we 
are entitled to take the EdS background without 'Dark 
Energy' and interpret the negative averaged curvature 
contribution as a replacement of the 'Dark Energy' in 
the standard interpretation, where curvature is strictly 
zero on all scales. We have quantified this contribution 
as less than 20% for one-sigma fluctuations on 200 Mpc, 
which is - in this model - not enough to compensate the 
missing 'Dark Energy' of the order of 73%. (Note that 
when speaking about 'Dark Energy' we do not imply that 
the X— Matter component always produces an accelerat- 
ing universe model.) For two-sigma fluctuations on this 
scale we find 35% of the needed amount, but due to the 
architecture of our model, as discussed below, this con- 



tribution falls off rapidly on very large scales. We may 
compare this regional behavior in the sense of so-called 
'void models' that assume we are living in a large under- 
dense region, see e.g. ]7, 8, 43] and references therein. 
In the spirit of these models our generic model 'explains 
away' 1/2 of the 'Dark Energy' for a slightly untypical 
region (two-sigma) on the scale of 200 Mpc. 

Before we explain why we cannot expect this behav- 
ior on larger scales in our model, we briefly sketch the 
opposite behavior of the backreaction on smaller scales, 
where a strongly anisotropic collapse into sheet-like and 
filamentary-like structures features a negative backreac- 
tion term Qd and a positive curvature term (TV) -p. This 
situation occurs e.g. on the scale of 50 Mpc, below which 
the X— Matter component starts to compete with the 
matter density and mimics the presence of 'Dark Matter'. 
This result points to the high relevance of the backreac- 
tion effect for the explanation of 'Dark Matter' in over- 
densities, while this effect decays on large scales leaving 
only a few percent remnant on an assumed homogeneity 
scale of about 400 Mpc. This latter is comparable with 
the 'Dark Energy' remnant on such scales. Both effects 
thus are small but still alter the large-scale cosmological 
parameters at the percent level, see Fig. 8. 



B. Conclusions 

There are several serious shortcomings of our model 
that has a limited architecture as compared with the gen- 
eral situation, below discussed as the background problem, 
and also our model is, due to the choice of CDM and 
ACDM initial conditions, exploited in a regime where 
the backreaction effect is not yet effective on large scales, 
below discussed as the amplitude problem. 



1. The background problem 

In our model the background is "fictitious": while our 
backreaction models have a generic structure, our quan- 
titative interpretations largely depend on the choice of 
background [40]. The reason is that both of our back- 
grounds contain a large amount of 'Dark Matter', while 
the effect we study produces kinematical contributions 
that act in a similar way. The same argument holds for 
the background containing 'Dark Energy'. The imple- 
mentation of a physical background, defined through the 
spatial average of the considered fluctuations, is the next 
step to be envisaged (compare the flrst results obtained 
in [58]). Such a physical background interacts with the 
backreaction models and thus changes their interpreta- 
tion. While in the corresponding Newtonian model [BKS] 
the background forms the average distribution due to its 
torus architecture ]10], i.e. the scale factor = a on the 
periodicity scale, the relativistic model studied here must 
be considered a hybrid construction featuring a volume- 
scale factor ax> and a background scale factor a that only 
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approximately matches the volume-scale factor on the 
largest scales as a result of large-scale remnants of back- 
reaction. An indication for this shortcoming is also that 
the backreaction model, seen on the FLRW background, 
only covers very restricted subcases of exact GR solu- 
tions, unlike the situation in Newtonian cosmology, where 
the plane and spherical collapse is in general covered by 
the backreaction model. 



2. The amplitude problem 

Cloesely related to the first problem of a background 
containing a dominating 'Dark Matter' component, we 
here point out a serious mismatch between the time of 
onset of the actual large-scale backreaction effect and the 
age of the Universe as it is determined by the background 
model. We can take as an indicator the onset of an ac- 
celerating period that will set in as a result of dominance 
of backreaction effects over the density. Such an acceler- 
ation may not be needed to explain observational data; 
also, it is a strong requirement in view of the smallness 
of the kinematical backreaction term, which is the only 
backreaction component entering the volume acceleration 
equation, while the impact of X-matter is considerably 
stronger. From our model we infer that this acceleration 
phase is seen in the future on scales roughly above 40 
Mpc. This result confirms the amplitude problem already 
reported in [64] , if we suppose that backreaction were to 
replace 'Dark Energy' completely. Also, if we would con- 
sider the analyses presented in this paper in the period 
of onset of the effect, our quantitative conclusions on the 
amount of AT— Matter would substantially change. Thus, 
all of our results remain conservative - although visible 
in the fiuctuation properties on intermediate and small 
scales - since the actual effect is postponed to the future 
due to our standard initial data setting. 

We can understand this by a simple consideration: the 
averaged density decays as ^ the constant curvature 
(contained in the full averaged Ricci curvature) decays 
as o^'j^ , 

the leading large-scale mode of the backreac- 
tion model (the kinematical backreaction Qp) decays as 
and, finally, the non-decaying cosmological constant. 
Thus, we expect several periods of onset of dominance of 
(i) the constant-curvature part and later (ii) the backre- 
action contribution, assuming that the cosmological con- 
stant is set to zero. As we saw, this latter period lies in 
the future for our initial data setting. 

There are furthermore limitations due to the model's 
input of a first-order deformation. Obviously, this lim- 
its the range of applicability for the full backreaction 
regime that here lies in the future, but would anyway 
not be covered by our extrapolation of a first-order so- 
lution. Generally, our model is applicable only to large 
scales, well beyond the scales of virialized objects. The 
matter distribution is assumed to be fairly smooth and 
does not take into account the small-scale discretization. 
The importance of the modeling of small-scale struc- 



ture has to be emphasized. For example, the determi- 
nation of the actual volume fraction in devoid regions is 
needed to precisely quantify the contribution of a nega- 
tive curvature (and in turn of the contribution to what 
would be interpreted as 'Dark Energy' in the standard 
model), [52. 64, 66[, and this sensibly depends on the 
high-resolution modeling of small-scale structure. Col- 
lapse models, as the one based on [BKS[, [38[; have to 
be refined by e.g. including velocity dispersion. A cor- 
responding argument, making use of the evolution of the 
fraction of virialized objects together with an attempt to 
also model the light-cone effect, may give an indication 
of its quantitative importance [56[ . 



3. Outlook 

Clearly, in view of the remarks above, further major 
efforts are needed, and the present investigation can only 
be considered as a further step towards quantifying back- 
reaction effects. We expect, however, that the elements 
we provided in this work will prove useful for further 
considerations and the construction of improved models. 
At any rate, a model based on the RZA is expected to 
realistically model the large-scale skeleton, and so lies 
at the heart of any model for structure formation includ- 
ing fully relativistic numerical simulations. Furthermore, 
the above remark has to be seen in junction with another 
major effort to reinterpret observational data in the rel- 
ativistic models: the subtleties and complexity of the in- 
terpretation of curvature energies as 'Dark Energy' and 
'Dark Matter', seen on intermediate scales in our inves- 
tigation, are not only a challenge in view of the scale- 
dependence of the effect, but imply important changes of 
light-propagation in these models [31, 44, 55[. Here again 
the precise modeling of small-scale structure is crucial. 

The present investigation reveals, in a quantified way, 
how we have to proceed in order to master inhomoge- 
neous universe models. We identify four key-questions 
to be addressed: first, how can we deal with the 'kine- 
matical Dark Matter' produced in the relativistic models 
on all scales, and how does this alter the initial conditions 
and their compatibility with CMB constraints? Second, 
how can we improve the model to quantify the effect on 
smaller scales than those accessible here? Third, how 
can we implement a perturbation scheme on the physical 
background, given by the average of the model? Fourth, 
how are observational data reinterpreted in the relativis- 
tic models? 

The answer to the first three questions would alter the 
time-scale of the background, of structure formation and 
of our interpretation of backreaction. Would this imply a 
shift of the backreaction regime into an epoch before to- 
day? Interpreting the onset of an acceleration period as 
coinciding with an apparent acceleration is highly prob- 
lematic [6[, we are tempted to require such a shift to 
solve the 'coincidence problem', but an acceleration of 
the model may after all not be needed to explain ob- 
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servational data. The answer to the third question will 
deliver backreaction models that interact with the back- 
ground; they will contain, e.g., the volume scale factor 
(and not the FLRW scale factor), which itself depends 
on the backreaction terms. An iterative scheme results 
that would substantially modify the time-scales of our re- 
sults. In turn, the amplitude of large-scale backreaction 
could substantially exceed the 6% effect already present 
in our model on an assumed homogeneity scale of 400 
Mpc. 

Finally, let us emphasize again that we constructed 
our model by keeping basic cornerstones of the standard 
model in the early stages of the cosmic evolution. In- 
sufficiency of late-time backreaction effects may indicate 
that a more drastic paradigmatic change is needed that 
takes into account backreaction effects in the Early Uni- 
verse: consequences of emerging average properties from 
inhomogeneous inflationary models [20], and cosmologi- 
cal models that are built on global principles like globally 
stationary cosmologies [14], which inherently explain a 
large backreaction component contained in the physical 
background, have thus far not been thoroughly analyzed. 
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7i = 72 = 73 = 0, what, by the definition of Eq. (50), is 
equivalent to the right-hand-side of Eq. (Al). □ 



Proposition 2. 



1 



(^■^^n(eg)^ = ±r^i(ev)); 
' (""^ni(eg)^^l(«-Ai(ev)); 

(A3) 

Proof: 

• Once again the vanishing of the backreaction term, 
given the right-hand-side, is straighforward. 

• For the other way we use (Prop. 1): if, Vt, ^^^Q© = 0, 
then one can write (Ili)j and (Illi)^; as functions of (li) j; 
with Eqs. (45) and (48) we have: 



(A4) 

where X = ^ /(I + - This finally leads to the 

3 

equalities of the right-hand-side of Eq. (A3). □ 



Appendix B: Orthonormal basis representation 



Appendix A: Demonstration of Propositions 1 and 2 
Proposition 1. 

r (iii)^ = 1 iwi , 

"^^^Qv = 0^1 (Al) 
[ (IIIi)^ = l(Ii)| . 

Proof: 

• With Eq. (50), the demonstration of the vanishing 
of the backreaction term, given the right-hand-side, is 
straighforward. 

• For given vanishing backreaction, i.e. supposing that: 
Vt, ^^^Qt> = we proceed as follows. Excluding the 
case ^(i) ~ Vt, using Eq. (50), the vanishing of the 
backreaction term is equivalent to: 

Vi, 7i+e(072+e'(073 = 0. (A2) 

As the monomials in ^ are linearly independent, 
already the coefficients have to vanish. Therefore, 



We will, in this appendix, give all the relevant expres- 
sions for the case where we employ the standard assump- 
tion of orthonormal frames (Option 1 in the text, see 
III A 2) . The coframes are in this case given by 

{t, X^-) := a {t) (j^J + P", + at)P") , (Bl) 

which combine, via Sab, to the complete metric coeffi- 
cients 

g^r-^ SabV^^v) ■ (B2) 

We evaluate for this case some relevant fields furnishing 
the RZA and its average properties: 

• the coefficients of the metric tensor: 

"^^^g^.it.X") = a\t) {% + 2P(y) + Pk^P'^j 

+e{t)Pk^P''j} ■ (B3) 
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• the local volume deformation: 



1 



6 



J := -eabce--'-TtV''jV''k 



a'm ■ (B4) 
• the first scalar invariant of the expansion tensor: 



(B5) 



• the second scalar invariant of the expansion tensor: 



\a{t)) ^a{t)Z 

2 1 a mz] ■ 



(B6) 



• the third scalar invariant of the expansion tensor: 
1 



a{t)j ^\a{t)) 3 

ia{t) fi ^{t)i\ , ^^{t)a 



(B7) 



ddd 



2aW \Z iit)Zj 3 
We have introduced the functions: 

X') := 5i + CAi + e'^dd + e'^ddd , 

Si -.^ S{ti, X'') := f + li + Hi + Illi , 



Dd:^Dditi,X'') := Id (1 + 1; + 11;) + P'.P^feP 
Dm ■■= D^itu X'') := Ild (1 + li) + P'jP\P\ 



ryk 



Dddd '■= Dddditi, X'') :— Illd , 



li HP'',) 
ni:=II(P'^) : 
III; := : 

Id 

lid n(P",) 

Illd := III(P°) 



2'^"''='^ P k ^ 



1 

6 



-f 1, f*^*^ p"^ pb . pc, 



1 



6 



t , e'-i^ P" P'' 



Using Eq. (45) one can express the scalar invariants of 
the peculiar-expansion tensor as a function of J and ^ : 



^ 6 Li 3i 2^(t) U CW^. 



(B9) 

Inserting Eqs. (B8) into Eq. (47), the backreaction term 
in the 'Relativistic Zel'dovich Approximation' reads: 



RZA 



(BIO) 



or finally, using the first equation of Eqs. (B8): 



RZAg 



•D 



(71 + ^72 + 



((^i)x + ^ (^d)x + e (^dd)x + ?3 (Dddd)x)' 



with: 



-2(^i)i(^dd>x-|(^d>x, 



6 J (i?ddd)x - I (^d>x (■Ddd)^ 

2 /T-. \2 



X \^ddd/x 3 \-^d/x \-^dd/x : 

= 2(i?d)x(^ddd>x- |(-Ddd>x- 



(Bll) 
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